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We present a search for pair production of doubly-charged Higgs bosons in the processes qq — > 
H ++ H decaying through H — > t ± t ± , fj, t , fi ± [i ± . The search is performed in pp collisions at 
a center-of-mass energy of y/s = 1.96 TeV using an integrated luminosity of up to 7.0 fb _1 collected 
by the DO experiment at the Fermilab Tevatron Collider. The results are used to set 95% C.L. limits 
on the pair production cross section of doubly-charged Higgs bosons and on their mass for different 
H branching fractions. Models predicting different H decays are investigated. Assuming 
B(H ±± —> t t ) = 1 yields an observed (expected) lower limit on the mass of a left-handed 
boson of 128 (116) GeV and assuming B{H — > fi ± r ± ) = 1 the corresponding limits are 144 
(149) GeV. In a model with B(i3' ±± -> r ± r ± ) = B(H ±:t -> (j^t*) = B(H ±± -> t l± t L± ) = 1/3, we 
obtain Af(i/± ± ) > 130 (138) GeV. 

PACS numbers: 14.80.Fd,13.85.Rm 



Doubly-charged Higgs bosons (iJ ±=t ) appear in models 
with an extended Higgs sector such as the Little Higgs 
model left- right symmetric models Q, and in mod- 
els with SU(S) c xSU(3) L xU(l) Y (3-3-1) gauge symme- 
try 1- 

The H ±=t bosons could be pair-produced and observed 
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at a hadron collider through the process qq — > Z/j* — > 
H++H— -> £ + £' + £-£'~ (£,£' = e,/*,r). Single pro- 
duction of 77 ±± bosons through W exchange, leading 
to i? ±± i? T final states, is not considered in this Let- 
ter to reduce the model dependency of the results 0]. 
Some models favor a mass of the i? ±± boson at the elec- 
troweak scale ■ The decay into like-charge lepton pairs 
violates lepton flavor number conservation. The decays 
H ±:t — > r ± r ± are predicted to dominate in some scenar- 
ios, such as the 3-3-1 model of Ref. [f|. In a Higgs triplet 
model that is based on a seesaw neutrino mass mech- 
anism, a normal hierarchy of neutrino masses leads to 
approximately equal branching fractions for H ±± boson 
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decays to tt, fj/r, and if the mass of the lightest neu- 
trino is less than 10 meV Q. In this Letter, we present 
the hrst comparison of data with this model and the first 
search for if ±± — > t ± t ± decays at a hadron collider. 

In left-right symmetric models, right-handed states 
(H^) appear in addition to left-handed states (H^). 
They are characterized through their coupling to right- 
handed and left-handed fermions, respectively. The cross 
section for production of right-handed H 
about a factor of 2 smaller than for H^ + H~ 
of the different coupling to the Z boson 



pairs is 
L because 
. The mass 



limits for bosons therefore tend to be weaker than 

for bosons. 

Searches for production of H bosons have been per- 
formed previously at the CERN e+e" Collider (LEP) @ 
and at the DESY ep Collider (HERA) [jjj. Limits on 
the mass of the iJ ±± boson were obtained in the range 
of 95 — 100 GeV, depending on the flavor of the final state 
leptons. The OPAL and HI Collaborations searched 
for single H ±zt production in the processes e + e~ — > 
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and e ± p 



and through 
con- 



11] 



the study of Bhabha scattering e + e~ — » 
straining the H ±±: boson's Yukawa couplings h ee to elec- 
trons. Bounds on decays such as r — > or /i — > e^f and 
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the measured (g — 2)^ also constrain different h 
At the Fcrmilab Tevatron Collider, the DO and CDF Col- 
laborations published limits for /i^t, ee, er, and fir final 
states in the range M(H^) > 112 — 150 GeV, assuming 
100% decays into the specified final state [l3l - fl6j . 

The results in this Letter are based on data collected 
with the DO detector at the Fermilab Tevatron Collider 
and correspond to an integrated luminosity of up to 
7.0 fb _1 . The DO detector [17J comprises tracking detec- 
tors and calorimeters. Silicon microstrip detectors and a 
scintillating fiber tracker are used to reconstruct charged 
particle tracks within a 2 T solenoid. The uranium and 
liquid-argon calorimeters used to measure particle ener- 
gies consist of electromagnetic (EM) and hadronic sec- 
tions. Muons are identified by combining tracks in the 
central tracker with patterns of hits in the muon spec- 
trometer. Events are required to pass triggers that select 
at least one muon candidate. 

All background processes are simulated using Monte 
Carlo (MC) event generators, except the multijet back- 
ground, which is determined from data. The T^+jet, 
Z/-f* — > £ + £~, and tt processes are generated using 
ALPGEN [18| with showering and hadronization provided 
by pythia [l9| . Diboson production ( WW, WZ, and ZZ) 
and signal events are simulated using pythia. The signal 
samples for the model with equal branching ratios for the 
decays H ±:t — > t ± t ± , /M^fJ^, and ^ ± r ± are generated 
using Yukawa couplings h^ T = h Tfl = \/2h TT = \/~2h 
The tau lepton decays are simulated with tauola [2 
which includes a full treatment of the tau polarization 
All MC samples are processed through a GEANT 2l| sim- 
ulation of the detector. Data from random beam cross- 



ings are overlaid on MC events to account for detector 
noise and additional pp interactions. The simulated dis- 
tributions are corrected for the dependence of the trigger 
efficiency in data on the instantaneous luminosity and 
for differences between data and simulation in the recon- 
struction efficiencies and in the distribution of the lon- 
gitudinal coordinate of the interaction point along the 
beam direction. Next-to-leading order (NLO) quantum 
chromodynamics calculations of cross sections are used 
to normalize the signal and the background contribution 
of diboson processes, and next-to-NLO calculations are 
used for all other processes. 

Two types of tau lepton decays into hadrons (r^) are 
identified by their signatures: Type-1 tau candidates con- 
sist of a calorimeter cluster, with one associated track 
and no subcluster in the EM section of the calorimeter. 
This signature corresponds mainly to r ir^v decays. 
For typc-2 tau candidates, an energy deposit in the EM 
calorimeter is required in addition to the typc-1 signa- 
ture, as expected for r ± — > ^tx^v decays. The outputs 
of neural networks, one for each tau- type, designed to 
discriminate t% from jets, have to be NN T > 0.75 [22l |. 
Their input variables are based on isolation variables for 
objects and on the spatial distribution of showers. The 
tau lepton energy is measured with the calorimeter. 

We select events with at least one muon and at least 
two Tfr candidates. The muons must be isolated, both 
in the tracking detectors and in the calorimeters. Each 
event must have a reconstructed pp interaction vertex 
with a longitudinal component located within 60 cm of 
the nominal center of the detector. The longitudinal 
coordinate Zdca of the distance of closest approach for 
each track is measured with respect to the nominal cen- 
ter of the detector. The differences between z^ ca of the 
highest-pr muon and the two highest-pr th (labeled T\ 
and r 2 ), must be less than 2 cm. The pseudorapidity Q 
of the selected muons, n, and t 2 must be |f? M | < 1.6 
and |r; T1 ' 2 | < 1.5, respectively, and for additional Th can- 
didates we require |t7 t | < 2. The transverse momenta 
must be p% > 15 GeV and p^ 1 ' 2 > 12.5 GeV. All se- 
lected Th candidates and muons arc required to be sep- 
arated by AK^r > 0.5, where AK = (A</>) 2 + (Arj) 2 
and (j) is the azimuthal angle, and the two leading Th 
must be separated by AlZ TlT2 > 0.7. The sum of the 
charges of the highest-pr muon, n , and r 2 is required to 



be Q = Zi 



±1 as expected for signal. After 



all selections, the main background is from diboson pro- 
duction and Z — > r + r _ , where an additional jet mimics 
a lepton. 

We estimate the multijet background using three in- 
dependent data samples and identical selections, except 
with the NN T requirements reversed, by requiring that 
either one or both Th candidates have NN T < 0.75. The 
simulated background is subtracted before the samples 
are used to determine the differential distributions and 
normalization of the multijet background in the signal 
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FIG. 1: (color online). A'I(ti,T2) distribution for the (a) q T1 = q T2 and (b) q T 



-q T2 samples, and (c) transverse momentum 



of the doubly-charged dilepton system p T , for all four samples combined, after all selections. The data are compared to the 



sum of the expected background and to simulations of a signal for M{H ) 



= 1, and B{H^ 



=) = B(H±± 



M*) = B{Hf- 



120 GeV and B(H^ 



±± 



= 1, 



calculation of the cross section. "Other" background comprises W+jet, Z/'y* 
the range of the histogram are added to the last bin. 



) = 1/3, normalized using the NLO 
, and tt processes. All entries exceeding 



TABLE I: Numbers of events in data, predicted background, 
and expected signal for M(H^ ) = 120 GeV, assuming the 
NLO calculation of the signal cross section for B(H^ — > 
r ± r ± ) = 1, B(H^ -> ^T*) = 1, and B(H^ -> r ± r ± ) = 
BiH^ -» fj, ± fj, ± ) = B{H^ -»■ h ± t ± ) = 1/3. The numbers 
are shown for the four samples separately, together with their 
total uncertainties. 
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Background 
















Z -> T + T~ 


8.2 ± 1.1 


3.4 


± 


0.5 


4.8 ±0.7 


< 0.1 


< 0.1 


Z — > /i + £t~ 


5.1 ±0.7 


2.2 


± 


0.3 


2.5 ±0.4 


0.1 ±0.1 


0.2 ± 0.1 


Z -> e+e- 


0.3 ±0.1 


< 





1 


0.3 ±0.1 


< 0.1 


< 0.1 


W + jets 


2.9 ±0.4 


1.1 


± 


0.2 


1.8 ±0.3 


< 0.1 


< 0.1 


it 


0.6 ±0.1 


0.3 


± 


0.1 


0.3 ±0.1 


0.1 ±0.1 


< 0.1 


Diboson 


10.5 ± 1.7 


0.5 


± 


0.1 


8.5 ± 1.4 


0.4 ±0.1 


1.1 ± 0.2 


Multijet 


< 0.8 


< 





2 


< 0.5 


< 0.1 


< 0.1 


Background 
















Sum 


27.6 ±4.9 


7.5 


± 


1.2 


18.2 ±3.3 


0.6 ±0.1 


1.3 ± 0.2 


Data 
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region. A second method used to estimate the multijet 
background is based on the fact that events with Q = ±1 
are signal-like, whereas events with Q = ±3 correspond 
largely to multijet background. To reduce the VF+jets 
contribution in the sample with Q = ±3, the visible 
W boson mass Mw = \J^V^"4t{X ~ cos</>) is required 
to be < 50 GeV, where p^ is the muon momentum, $t 
the imbalance in transverse momentum measured in the 
calorimeter, and <fi is the azimuthal angle between the 
muon and the direction of the $t- The total rate of ex- 



pected multijet background events following all selections 
is negligible (< 3% of the total background). We also use 
the sample where both 77, candidates have NN T < 0.75 
to study the rate of jets that arc falsely reconstructed 
as Th and we find this rate to be well modeled by the 
simulation. 

To improve the discrimination of signal from back- 
ground, the data are subdivided into four nonoverlapping 
samples, depending on the charges of the muon and 
the Th candidates (q T ) and the number of muons (N^) 
and Th (N T ) in the event. First, we define two samples 
for events with = 1 and N T = 2. Because the two 
like-charge leptons are assumed to originate from a sin- 
gle if ±=t decay, we consider separately events where both 
tau leptons have the same charge, q Tl = q T2 , and events 
with T\ and T2 of opposite charge, i.e., q Tl = —q T21 which 
implies that one of the r leptons and the muon have the 
same charge. The third sample is defined by N T = 3 and 
the fourth sample by = 2, without any additional 
requirements on the charges. 

The distributions of the invariant mass of the two lead- 
ing tau candidates, M(ti,T2), for the like and opposite- 
charge samples are shown in Figs. [IJa) and (b). The 
separation into samples with different fractions of sig- 
nal and background events increases the sensitivity to 
signal, as the composition of the background is differ- 
ent, with the like-charge sample being dominated by 
background from Z+jets decays and the opposite-charge 
sample by background from diboson production. The 
diboson background is mainly due to WZ — > [ii>e + e~ 
events where the electrons are misidentified as tau lep- 
tons. In Fig. [TJc) we show the transverse momentum 
of the doubly-charged dilepton system, pj,, which corre- 
sponds to the reconstructed if ±=t — s- £ ± £' ± decay, where 
£ ± £' ± = (//^Tjk, /i ± T 2 ± , T-j^r^) is the pairing of the two 
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FIG. 2: (color online). Upper limit on the H L 
production cross section for (a) B(H^ 
(b) B{H^ -> M ±r± ) = 1, and (c) B(H^ -+ r ± r ± ) = 
BiH^ -> fit fi*) = B{H^ -> /i ± r ± ) = 1/3. The bands 
around the median expected limits correspond to regions of 
±1 and ±2 standard deviation (s.d.), and the band around 
the predicted NLO cross section for signal corresponds to a 
theoretical uncertainty of ±10%. 



highest-pT T~h and the highest-pT muon that have the 
same charges. Since \Q\ = 1, only one such pairing exists 
per event. The expected number of background and sig- 
nal events for the four samples and the observed numbers 
of events in data are shown in Table |T] with the statistical 
uncertainties of the MC samples and systematic uncer- 
tainties added in quadrature. 

Since the data are well described by the background 
expectation, we determine limits on the H ++ H pro- 



TABLE II: Expected and observed limits on M(H ) (in 
GeV) for left-handed and right-handed H bosons. Only 
left-handed states exist in the model that assumes equality of 
branching fractions into rr, fir, and /i/i final states. We only 
derive limits if the expected limit on M(H ±± ) is > 90 GeV. 



Decay 
B(H ±:k -> T+T*) = 1 

B(H ±:t -> i^t*) = 1 
Equal B into 
T t , /i ± ^i ± , fj, T 




expected observed expected observed 
116 128 

149 144 119 113 
130 138 


-> fJ^/J?) = 1 


180 168 154 145 



duction cross section using a modified frequentist ap- 
proach 24J. A log-likclihood ratio (LLR) test statistic 
is formed using the Poisson probabilities for estimated 
background yields, the signal acceptance, and the ob- 
served number of events for different ff ±:t mass hypothe- 
ses. The confidence levels are derived by integrating 
the LLR distribution in pseudoexperiments using both 
the signal-plus-background (CL s+ b) and the background- 
only hypotheses (CLb). The excluded production cross 
section is taken to be the cross section for which the con- 
fidence level for signal, CL S =CL s+ i,/ CL},, equals 0.05. 
The M(ji,T2) distribution is used to discriminate signal 
from background. 

Systematic uncertainties on both background and sig- 
nal, including their correlations, are taken into account. 
The theoretical uncertainty on background cross sections 
for Z/Y* £ + £~, W+jets, tt, and dibo son production 
vary between 6% — 10%. The uncertainty on the mea- 
sured integrated luminosity is 6.1% 25[. The systematic 
uncertainty on muon identification is 2.9% per muon and 
the uncertainty on the identification of Th, including the 
uncertainty from applying a neural network to discrimi- 
nate Th from jets, is 4% for each type-1 and 7% for each 
type-2 Th candidate. The trigger efficiency has a system- 
atic uncertainty of 5%. The uncertainty on the signal 
acceptance from parton distribution functions is 4%. 

In Fig. [2J the upper limits on the cross sections are 
compared to the NLO signal cross sections for 
pair production 0] for some of the branching ratios con- 
sidered. The corresponding expected and observed limits 
are shown in Table HH 

The boson mass limits assuming B(H^ — > 

r ± r ± ) + B(H ±:t —> /i ± /i ± ) = 1 are determined by com- 
bining signal samples generated with pure 4r, (2t/2/z), 
and 4/x final states with fractions £> 2 , 223(1 — B), and 



(1 — B) , respectively, where B = B(H 



±± 



Here, we include in the limit setting the distribution of 
the invariant mass of the two highest px muons, includ- 
ing the systematic uncertainties and their correlations, 
from a search for H ++ H — > 4/i decays performed by 
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(b) DO, L up to 7.0 fb 1 
pp^ HfHf 
B 8S = B„. = 
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NKH**) (GeV) 

FIG. 3: (color online). Expected and observed exclusion re- 



gion at the 95% C.L. in the plane of B(H 



7 7 ) versus 



M(H ±:t ), assuming B(H ±:t -> r ± r ± )+B( J ff ±± 
1, for (a) left-handed and (b) right-handed i/ ±=t bosons. The 
band around the expected limit represents the uncertainty on 
the NLO calculation of the cross section for signal. 



the DO Collaboration in 1.1 fb 1 of integrated luminos- 
ity (l3j . The results are shown in Fig. [3] for varying 
B = 0%- 100% in steps of 10%. When performing this 
analysis, we found that the statistical uncertainties on 
the background simulations were overestimated in (]~3j ) . 
A standard treatment of the uncertainties in the limit 
setting improves the mass limits for the 4[i final state, as 
shown in Table [TTJ 

In summary, we have performed the first search at 
a hadron collider for pair production of doubly-charged 
Higgs bosons decaying exclusively into tau leptons. We 
set an observed (expected) lower limit of M(H L ) > 



±± 



128 (116) GeV for a 100% branching fraction of H 
r ± r ± , M(flf =*=) > 144 (149) GeV for a 100% branching 
fraction into /Ltr, and M( J ff^ ± ) > 130 (138) GeV for a 
model with equal branching ratios into rr, /it, and mi. 
These are the most stringent limits on H boson masses 
in these decay channels. 

We thank the staffs at Fermilab and collaborating 
institutions, and acknowledge support from the DOE 
and NSF (USA); CEA and CNRS/IN2P3 (France); 
FASI, Rosatom and RFBR (Russia); CNPq, FAPERJ, 
FAPESP and FUNDUNESP (Brazil); DAE and DST (In- 



dia); Colciencias (Colombia); CONACyT (Mexico); KRF 
and KOSEF (Korea); CONICET and UBACyT (Ar- 
gentina); FOM (The Netherlands); STFC and the Royal 
Society (United Kingdom); MSMT and GACR (Czech 
Republic); CRC Program and NSERC (Canada); BMBF 
and DFG (Germany); SFI (Ireland); The Swedish Re- 
search Council (Sweden); and CAS and CNSF (China). 



[1 
[2 

[3 

[5, 

[6 
[' 

[8 
[9 

[10 

[11 
[12 

[13 
[14 
[15 
[16 

li- 
tis 

[19 
[20 
[21 



N. Arkani-Hamed et al, J. High Energy Phys. 08, 021 

(2002) . 

R.N. Mohapatra and G. Senjanovic, Phys. Rev. Lett. 44, 
912 (1980); R.N. Mohapatra and G. Senjanovic, Phys. 
Rev. D 23, 165 (1981); J.F. Gunion et aJ., Phys. Rev. 
D 40, 1546 (1989); N.G. Deshpande et al, Phys. Rev. D 
44, 837 (1991). 

J.E. Cieza Montalvo et al, Nucl. Phys. B756, 1 (2006); 
erratum- ibid. B796, 422 (2008). 

A.G. Akeroyd and M. Aoki, Phys. Rev. D 72, 035011 
(2005). 

C. S. Aulakh, A. Melfo, and G. Senjanovic, Phys. Rev. D 
57, 4174 (1998); Z. Chacko and R. N. Mohapatra, Phys. 
Rev. D 58, 015003 (1998). 

E. Ramirez Barreto, Y.A. Coutinho, and J. Sa Borges, 
Phys. Rev. D 83, 075001 (2011). 

M. Kadastik, M. Raidal, and L. Rebane, Phys. Rev. D 
77, 115023 (2008); A. Hektor et al, Nucl. Phys. B787, 
198 (2007). 

M. Miihlleitner and M. Spira, Phys. Rev. D 68, 117701 

(2003) and private communications. 

G. Abbiendi et al. (OPAL Collaboration), Phys. Lett. 
B 526, 221 (2002); P. Achard et al. (L3 Collaboration), 
Phys. Lett. B 576, 18 (2003); J. Abdallah et al (DELPHI 
Collaboration), Phys. Lett. B 552, 127 (2003). 
A. Aktas et al (HI Collaboration), Phys. Lett. B 638, 
432 (2006). 

G. Abbiendi et al. (OPAL Collaboration), Phys. Lett. B 
577, 93 (2003). 

J.F. Gunion et al, Phys. Rev. D 40, 1546 (1989); 
R.N. Mohapatra, Phys. Rev. D 46, 2990 (1992). 
V.M. Abazov et al. (DO Collaboration), Phys. Rev. Lett. 
101, 071803 (2008). 

V.M. Abazov et al. (DO Collaboration), Phys. Rev. Lett. 
93 141801 (2004). 

D. Acosta et al. (CDF Collaboration), Phys. Rev. Lett. 
93, 221802 (2004). 

T. Aaltonen et al. (CDF Collaboration), Phys. Rev. Lett. 
101, 121801 (2008). 

V.M. Abazov et al. (DO Collaboration), Nucl. Instrum. 
Methods Phys. Res. A 565, 463 (2006); M. Abolins et al., 
Nucl. Instrum. Methods in Phys. Res. A 584, 75 (2008); 
R. Angstadt et al, Nucl. Instrum. Methods in Phys. Res. 
A 622, 298 (2010). 

M.L. Mangano et al, J. High Energy Phys. 07, 1 (2003); 
we use version 2.11. 

T. Sjostrand, S. Mrenna, and P. Skands, J. High Energy 
Phys. 05, 026 (2006); we use version 6.323. 
Z. Was,, Nucl. Phys. B Proc. Suppl. 98, 96 (2001); we 
use version 2.5.04. 

R. Brun and F. Carminati, CERN Program Library Long 



8 



Writeup W5013, 1993. beam direction. 

[22] V.M. Abazov et al. (DO Collaboration), Phys. Rev. D 71, [24] W. Fisher, FERMILAB-TM-2386-E (2006). 

072004 (2005); erratum-ibid. D 77, 039901 (2008). [25] T. Andeen et al, FERMILAB-TM-2365 (2007). 

[23] The pseudorapidity is defined as r\ — — ln[tan(#/2)], 

where 6 is the polar angle with respect to the proton 



